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Summary 


BULLÖCK, D, A, (1975)—The general water circulation of Spencer Gulf, South Australia, in 
the period February to May. Trans. R. Soc. S. Aust. 99(1), 43-53, 28 February, 1975. 


Historical temperature and salinity data are presented together with the three-dimensional 
velocity structure of the general water circulation of Spencer Gulf, deduced from the data 
using a nUmerical model. In the southern area of the Gulf a cyclonic gyre is found which 
exchanges Gulf water wilh the ocean water outside the Gulf. On the western coast of the 
Gulf, a moderate boundary flow, which is called the “Port Lincoln Boundary Current”, is 
evident at all depths. The water in the northern portion of the Gulf circulates both in the 
vertical and the horizontal and appears to form an almost closed system separate from the 


system in the southern region, 


Introduction 

Ao investigation of the general water circu- 
lation in Spencer Gulf, and the region outside 
it, is made below from observation of property 
distributions and using a numerical model. The 
term “general water circulation” is used so that 
tidal currents which basically involve a periodic 
north-soulh movement of the Gulf’s waters are 
exeluded. Tidal currents can be superimposed 
on the more sustained flows caused by thermo- 
haline (and wind) forcing which are particu- 
larly important because of the part they play 
in the interchange of water in a substantially 
enclosed body of water. Thermohaline currents 
are brought about by horizontal pressure gra- 
dients due to density variations in the Gulf. 
which are caused by temperature and salinity 
differences, produced in part by outside water 
entering the Gull, and in part by variation in 
the effect of evaporation over the Gulf. 

All available data for the region, obtained 
on C,S.LR.Q. cruises, are shown in Fig. 1, The 
crises took place between 1961 and 19466 dur- 
ing the months indicated in the caption. No 
direct curren, readings are available, although 
some results on surface drift do exist from a 
drift card experiment (Marshallsay & Radok 
1972). 


Presentation and Interpretation of Cruise Data 

Using measurements made by F.R,V. Inves- 
figator (C.S.LR.O. Aust., 1968c) on 7 May 
1964, a vertical salinity section (Fig. 2a) was 
drawn for a line across the mouth of Spencer 
Gulf from Port Lincoln on the western side to 
Corny Point on Yorke Peninsula. The salinity 
is approximately homogeneously distributed 
over depth for any particular water colunm, 
although around the centre the situation is 
closer to that of a two-layer system since the 
salinity is almost constant with depth until 
about 20 m, at which point a steep salinity 
gradient occurs. Below this. the salinity is 
ugain constant with depth but greater than 
before, Also, the salinity increases with dis- 
tance fram west to cast, Indicating that water 
enters the Gulf on the western side and leaves 
via the eastern side after ils salinity has in- 
creased due to evaporation and mixing. The 
density section (Fig. 2b) has a similar form to 
the salinity section, indicating that salinity ts 
the dominant factor in determining thermo- 
haline flow. 

‘The most important feature of the hydrology 
of Spencer Gulf. evident from all data (Figs 
10, 17), is that the water is approximately 
homogencous fram top to bottom, and hence 
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Fig. 2. Vertical sections across Spencer Gulf near 
the mouth, Marked — - on Figure 1. 
data taken from CSIRO Aust. 1963, samp- 
ling depths indicates by dots, 
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it is reasonable that the quantities be repre- 
sented by their depfh-average, Fig, 3a—c shows 
horizontal depth-average distributions of salin- 
uy, temperature and density drawn from data 
obtained in the period 4-7 May 1964 
{C.S1,R.0. Aust, 1968c). Although other mea- 
surements have been made in Spencer Gulf, 
this cruise was by far the most comprehensive. 
From experience it was found necessary 1o 
restrict the data used to that obtained on one 
cruise, whieh spans only a short duration, so 
that seasonal and shorter term variations in 
properties do nol cause misleading interpreta- 
lions to be made. 


Fig, |. 
Key to Symbols 
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The distributions shown in Fig. 3 are con- 
sistent with those in Fig, 2 and give additional 
information. A flow up the western shore is 
evident and it can be seen that the salinity im- 
creases tastward over the whole southern area 
of the Guif indicating that the movement of 
water is basically clockwise. It appears that the 
adveclive exchange with the water outside the 
Gulf is limited to the area below a latitude of 
about 33°45’8. Above this region where the 
Gulf narows, the water appears “blocked-in”. 
Along the eastern and central areas of the 
upper part of the Gulf, the high salinity water 
seems to flow down the Gulf until it meets the 
lower salinity water taking part in the advec- 
tive process described above. 

The data from the northern part indicate that 
the water moves as a mass in phase with the 
tidal currents. At 34°00’S, 137° 10/B, two sta- 
lions (76 and 89) were occupied, separated by 
# period of 29 hours during which the depth- 
average salinity changed from 38.40%, to 
38.14%,. At 33°57'S, (37°25°E (stations 77 
und 85), the salinitites changed from 38.35%, 
to 38.60%; over a period of 21 hours. Using 
the horizontal salinity gradients at this location, 
Fig. 3 shows that this salinity variation corres- 
ponds to a tid excursion of a few kilometres 

The mechanism responsible for the homo. 
geneity of the water columns is based on ver- 
tical mixing promoted by tidal currents, which 
are amplified by the shallow water. running 
over a rough bottom which enhances vertical 
diffusion by the shear effect {Bowden 1965). 

The temperature variation over the Gulf is 
small (Fig, 3b); a plausible explanation for this 
observation is the following, Since, in the north 
of the Gulf, the advective exchange of pro- 
perlies is small and the water shallow, the 
effects of evaporation are more marked than 
further south, as is reflected by the high salini 


Location of Oceanographic Stations in the South Australian Gulf System until 1973, 


1962 (CSIRO Aust. 1967c} 
1967d) 


1966 (CSIRO Aust. 1969) 


19684) 


a HMAS. Gascoyne (G2/61) Feb.-March 1961 (CSIRO Aust 1966) 
© HMAS. Gasceyne (G2/63) March 1963 (CSIRO Aust. 19674) 
= ALM.A.S. Gascoyne (5/65) March-April 1965 (CSIRO 1967h) 
= HM.A.S, Gascoyne (G2/62) June-Aug. 

r HMAS. Gascoyne (G2/64) Feb. 1964 (CSIRO Aust, 

* HM_AS. Gascoyne (G2/65) Feb, 1965 (CSIRO Aust. 1968a) 
+ HMAS, Diamantina (Dm2/66) April 

2 ERV. Investigator 1963 (CSIRO Aust. (968b) 

* ERV. Investigator 1964 (CSIRO Aust. 1968c) 

< FRV. Investigator Dec. 1962 (CSIRO Aust. 

e R.V, Investigator Feb. 1965 (CSIRO Aust, 1968¢) 

= EV, Estelle Star April-May 1965 (CSIRO Aust. 1968f) 

» FY, Estelle Star April 1966 (CSIRO Aust. (968g) 

a ARV. Weerutta dan.-March 1961 (CSIRO Nake 1968h) 
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Vig. 3, Horizontal depth-average distributions in Spencer Gulf. Data taken from CSIRO Aust. 1968c. 
(a) Salinity (4%), (b) Temperature (°C), (e) s. 


ties. The process of evaporation requires the 
extraction of latent heat of vaporisation which 
in turn lowers the temperature and nullifies the 
elfect of direct heating. The following rough 
calculation shows that there is an approximate 
balance between the heat input from the atmos- 
phere and the heat output due to evaporation. 
From Haney (1972), the heat flux/unit area 
is proportional to the difference between the 
upparent air temperature and the sca tempera- 
ture (denoted by T,* and T,, respectively). 
Assuming this flux is alt available for evapora- 
tion, we obtain 


E= nA T) 


where E is the evaporative heat flux/unit area, 
and approximately à = 34 Wm=(°C)-! (Haney 
1972). Taking the values, 'T,* <= 21°C (the 
mean temperature for carly April), and E = 
100 Wnr?, and substituting in equation (1) 
yields T, =~ 18°C. which is approximately the 


observéd value (Fig. 3b), The neglect of advec- 
tive exchange is justifiable in view of the ob- 
served small variation of temperature and 
closed nature of the circulation in the northern 
Gulf, Vertical mixing maintains the tempera- 
ture difference between the atmosphere und the 
sea surface. 

Locations of the stations aré shown by the 
dots in Fig, 3a—c; stations were not occupied 
in some arcas, notably above Hardwick Bay on 
the eastern side of the Gulf. Nevertheless, the 
overall pattern of water properties can be in- 
ferred reasonably accurutely. 


The Numerical Model 


In order to investigate whether the above cir- 
culation could be considered as primarily due 
to thermohaline forcing or whether some other 
cause such as wind is dominant, and discover 
in greater detail the dynamics of Spencer Gulf, 
a numerical model was used to derive theore- 
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tically the field of flow from the density field 
shown in Fig, 3c. The model entitled FLOW- 
DEN computes the currents generated in a sea 
of density independent of depth, by wind and 
variations in atmospheric pressure, and by the 
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horizontal pressure gradients induced by the 
density feldt. This is done by solving the finite 
difference representation of the component 
transport equations written in a form applicable 
to the terrestrial situation”. 


au an H? dor , Tsx _ 2 
ne AHR eS RG, K| up] up (2) 
ae an H? do¢ , Tsy 

w gy y Pp akv sa 


where the following rotation is used; 


—H 


Ki 
udz, V= 
—H 


vdz. 


u = eastward component of velocity (east = x direction) 
v = northward component of velocity (north = y direction) 
up Ye = X & y components of bottom velocity (ug) 

z — distance measured perpendicular to earth's surface, positive upwards 
7, = (p—1)a 10°. pin kel 

f = Coriolis parameter = 2~cos?, w — angular speed of the earth’s rotation. # = latitude 
H ~ depth 
gê = [0# ey, nog, 235.00 = specific volume of the standard ocean 

g — gravitational acceleration 

n sca sutface elevation 

K = bottom drag coefficient 

mw — density of a standard ocean = — - 


The terms in equations (2) and (3) have the 
following physical meanings: the left-hand 
sides arc the linearized accelerations, the first 
terms on the right-hand sides are the Coriolis 
accelerations, the second terms are the bori- 
zontal pressure gradients due tọ elevation of 
the surface, the third terms are the thermo- 
haline terms Which take into account the hori- 
zontal pressure gradients due to the inhomo- 
geneity of the horizontal salinity and tempera- 
ture distributions. The final two terms are stress 
terms describing wind and bottom friction res- 
pectively- 

The magnitudes of the thermohaline terms 
were evaluated using another computer pro- 
gramme from the observed density distribution 
(Fig. 3c) and the depths were interpolated 
from those given on a naval hydrographic 
chart. A typical value of the x component of 


x & y components of surface wind stress 


the thermohaline terms was 300 mms. This 
can be compared with the typical magnitude of 
100 mms for the forcing term corresponding 
io wind stress, so thal the circulation in Spencer 
Gulf due to thermohaline causes can be cx- 
pected to be quite significant. The bottem drag 
coetiicient K = .0025. 

FLOWDEN was run with the precalculated 
thermohaline term as the only forcing term and 
with the thermohaline term together with a 
constant wind. After the equivalent of 17 hours 
an equilibrium current was reached (actually 
a small oscillation about the solution remained, 
but this was ignored). The initial conditions 
were that the How velocity was zero everywhere 
and the surface elevation was also zero every- 
where except at the boundary of the model 
actoss the south of the Gulf, which was 
assumed to have an initial slope similar to that 


1 Bye, J. A. T, 1975. Thallasso-scale numerical modelling. Computer Report. School of Earth Sciences. 
The Flinders University of South Australia. (Unpublished. ) 


2 Bullock, D: A. 1973, The general circulation of St. Vincent and Spencer Gulf, Honours Thesis, The 
Flinders University of South Australia. (Unpublished) 
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developed laler if the adjacent part of the 
Gulf. The boundary slope was found to have 
litle effect on the current distribution and so 
the assumption is presumed to nol have serious 
consequences, 

‘The resultant depth-aiverage current, which 
is the transport divided by the depth, was ob- 
tuned at each point of a grid system formed 
over the Gulf, in polar form so that the magni- 
tude and direction of the current could be seen 
immediately, The grid length was 6.4 km, 


rey) og UY BF Gy 
uz) = F+F Dy 
V * 30 


Applying equations (4) and (5) atz = a= 0 
and at z ~ —H gives the surface and bottom 
currents, 

Since », is independent of z it can be seen 
Irom equations (4) and (5) Ihat the vertical 
homogeneity of the Gulf implies that the cur- 
rent is a linear function of depth z. Thus a 
shear current exists. The right-hand side of 
equation (4) and (5) can be interpreted as fol- 
lows. The second terms are symmetric about a 
line marking mid-depth (Fig, 4) and give, as a 
function of depth, the current which is directly 
related to the focal distribution of density. The 
first terms on the right-hand side of the equa- 
tions, on the other hand, represent the depth- 
uverage current predicted by the model under 
the constraints imposed hy the houndaries of 
the Gulf (the shores and bed). 

The vertical components of velocity were 
also found, using the bottom current aml the 
gradients of the bottom topography, 


B OX B 3y 

Finally, the changes in the sea level which 
develop from the initial surface condition of 
zero slope are also obtained, 


i g (6) 


Results of the Numerical Mule 

(1) Predicted Velocity Fields due ta Thermo- 
haline forcing only 

Fig. 5 shows the main features of the pre- 
dicted depth-average current, The figure is con- 
sistent with the water mass analysis of the pre- 
vious section, ie. the circulation in the wide, 
southern part of Spencer Gulf consists ol a 
clockwise gyre, Water entering via the western 
side and leaving on the eastern side and the 
circulation in the nerthern end appears closed, 


To gain a detailed picture of the velocity 
field in the Gulf. a routine was added to the 
program which printed out the surface and bolt- 
tom currents. By involving the shallow waler 
assumption, i.e. that the Ekman layer (the sur- 
face layer influenced by wind) and the layer 
ol bollom frictional influence overlap (an 
assumption which is consistent with the waters 
vertical homogeneity), it can be shown (Bul- 
lock L973—footnole 2) thal the current vom- 
ponents at any level are given by, 


(+h C1) 
H 
(z + >) (w) 


The magnitude of the velocity yveeturs shown is 
typically about 0.07 ms? (ranging from vere 
ty O15 mst}, 


Prom Figs and 7, which show the surlace 
and bottom currents réspeciively, it can be scen 
thar at both the surface and bottum there is a 
very distinctive flow up the western side of the 
Gulf, which we shall call the “Port Lincoln 
Boundary Current”. This current is broader in 
extent ut the suufuce than at the bottom where 
it tends to peel away to the east down the slope 
of the bottom. The Port Lincoln Bounary Cur- 
rent is also somewhat stronger at ihe surfuce 
than at the hottom due to the frictional 
influence of the latter. 


The direction of the bottom current can he 
seen. to swing eastward and then ta the south- 
east over a broad arca before Mowing south- 
ward lo exit down the shovel-shaped topo- 
graphy through the Gull’s mouth. The surface 
current displays a similar pattern; however, 
since the Port Lincoln Boundary Current is 
wider al the surface, particularly near the 
mouth, the centre of the gyre of the surface 
circulation is located further to the north-eust 
than the centre of the’hottom gyre, Uvidence of 
the eastward movement of water from the 
boundary current is also provided hy Fig, 3u—c 
where the contours bulge to the east in this 
arei. 

The differences between the top and boltum 
current gyres implies that the surface and bot- 
tom currents on the eastern side of the Gulf 
are, in general, at an angle to cach other. un- 
like the western boundary. This explains why 
there is no obvious salinity or ¢ tongue indi- 
eating a return [Mow down the eastern side of 
the Gulf in Pig. 3b and c. 


WAIER CIRCULATION OF SPENCER GULF 49 


i} 

i 
A 
oa 
a 


aa 


$ 


Fig. 4. Current as a function of depth, with refer- 
ence to mid-depth. 


Another important feature to be noted from 
Figs 6 and 7 is that in the centre of the mouth 
the surface and bottom currents are in oppo- 
site directions so that a large shear is predicted 
here. The closeness of the density and salinity 
sections (Fig. 2b and 2a) to those of a two- 
layer system is also explained by this pheno- 
menon. It would be interesting to make direct 
current measurements in this area. 

Further north, where the Gulf starts to 
narrow (at about the latitude of Wallaroo}. the 
flow pattern is more complicated and conver- 


PORTAUGUSTA 


gences of the currents exist. A large proportion 
of the water in the northern end appears to 
recirculate in a vertical gyre. Fig. 8, which 
shows the vertical components of velocity, 
lends support to this intcrpretation by the 
downwelling at the top of the Gulf near Point 
Lowly, which corresponds to an inflexion in the 
salinity contours (Fig. 3a). In the southern 
part of the Gulf, a complex system ot vertical 
cells of typical dimension 20 km is predicted. 

The following changes in surface elevation 
are predicted by the model: in the southern 
part of the Gulf, a tow develops, the sea level 
on the western side being between 10 and 50 
mm higher than in the centre; the sea level also 
increases to the north. We can conclude that 
the circulation in Spencer Gulf which can be 
attributed solely to thermohaline causes is of 
significant magnitude and forms a major con- 
tribution to the water exchange processes of the 
Gulf, 


(2) Velocity field due to Thermohaline and 
Wind Forcing 

The numerical experiment was now extended 
to include the effect of wind stress in addition 


Fig. 5. Predicted depth-average thermohaline currents in Spencer Gulf. 
Fig. 6. Predicted surface thermohaline currents in Spencer Gulf. 
Fig. 7. Predicted bottom thermohaline currents in Spencer Gulf. 
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Fig. 8. Predieted regions of up- and down-welling due to thermohuline currents in Spencer Gulf. 
Fig. 9. Predicted currents in the presence of a 14 ms! SW wind. 


(a) Surface, (b) Bottom- 


ty the thermohaline forcing, The resulting sur- 
fuce and bottom currents predicted by the 
model for the same density strueture as before, 
combined with a 14 ms! south-west wind. are 
shown in Fig. 9a and 9b. 


The patiern of circulation is basically the 
sume but the velocity of the Now is accentuated, 
being typically about 0,16 ms, some velocity 
vectors having a magnitude aver 0.30 mst. 
Differences to the no-wind case do accur in the 
northern part of the Gulf where the surface 
eurrént is predominanily northward in direc- 
tion, particularly on the eastern side where i 
was directed southwards previously. Similarly 
the bottom current flows north adjacent to the 
shores, whereas in the no-wind case the bottem 
current was slight here. It can also be seen that 
with the wind added, the surface and bottom 
currents are almost the same everywhere. 


The model predicts that the surface elevation 
increasés greatly (up to 0,37 m) with distance 
up the Gulf, ie. the water is piled up at the 
northern end by the wind. The low in the 
southern half ts also accentuated, the Jevel 
being raised by up ta 80 mm around the shores 


compared with the depression of 10 mm in the 
centre, 


(3) Possible reasons for the Persistence of the 
Density Structure 

In the numerical model, the thermohaline 
forcing terms remained constant throughout the 
period fram zero velocity to equilibrium; Le 
we assume that the observed density field per- 
sists and is nat an instantaneous feature which 
subsequently runs down until the Gulf becomes 
a completely homogeneous body of water in 
which water motion of thermohaline origin 
would cease, 

It is surmised that the cause of the persist- 
ence of the thermohaline circulation is external, 
i.e. it is due to water from the Southern Ocean 
moving up under influences that determine ihe 
pattern of flow in this region, perhaps the pre- 
vuiling westerly winds. Fig. 12, which shows 
the dynamic topography of the ocean south of 
Australia (Bye 1971). indicates that the direc- 
tion of the geostrophic flow in the deep ocean 
is, in fact, basically perpendicular to the Aus- 
tralian coastline. 
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Fix. 10, Vertical sections along Investigator Strait 
ta the Continenial Slope. Marked —-— 
on Figure 1, data taken from CSIRO 
Aust. [966 and 1967a. Sampling depths 
indicated by dols. 

(a) Temperature (°C), (b) ee. 


The other factor which exerts a continuous 
effect on Spencer Gulf is evaporation, which 
has the effect of maintaining a high salinity and 
hence a high density at the northern end. It is 
worth noting that in Fig. 3a the salinity tongue 
extending downward from the Gulf’s northern 
end occurs over a favourable bottom gradient, 
Whereas this is not so for the tongue which 
corresponds ta the Port Lincoln Boundary Cur- 
rent. 


(4+) The Region Outside Spencer Gulf 

The locations of measurements made in the 
area outside Spencer Gulf are shown in Fig. 1. 
Data from this region were obtained primarily 
to aid the fishing industry and suffer from the 
drawback that there are no synoptic surveys of 
the whole area, which would enable horizontal 
distributions of properties to be drawn. Only 
vertical sections and profiles can be gatisfac- 
tarily constructed. The paths. along which sec- 
tions have been drawn are shown by the 
yarious lines in Fig. 1. 


$1 


Fig. 10a and b show temperature and +, 
sections drawn along a line extending through 
investigator Strait to past the end of the Con- 
tinental Slope, The datu were obtained in Feb,- 
March 1961 (C.S,LR,O. Aust. 1966) snd 
March 1963 (C.S.LR.0. Aust. 19674), One 
feature which is naticcable is the change in pint- 
tern from a vertically homogeneous situation 
just inside the mouth of Investigator Sirait to 
the more complicated horizontal stratification 
in the deeper waters outside the Gulf, 

The most significant feature shown by these 
sections, however, is the gradient in salinity, 
temperature and density about Station 122, 
which was located north-east of Cape Borda on 
the north-west tip of Kangaroo Island. This 
density gradient implies the existence of a 
thermohaline current, water to the west of Sta: 
tien 122 flowing approximately perpendicular 
to the plane of the section (into the page). This 
direction correlates very closely with the direc- 
tion of the thermohaline current entering 
Spencer Gulf which forms the Port Lincaln 
Boundary Current. Conversely there is also a 
southward directed flaw between Stations 121 
and 122 which could have its origins in the 
outhow from the eastern side of Spencer Gull. 
lt is possible that upwelling is also occurring 
near the longitude of Station 122. The peak in 
the density lines for the surface water under 
Station 41 is explained by the fact that this sta- 
tion was oecupied in a different year to the 
adjacent stations (although at the same time 
of the year), so that the surface water tempera. 
tures ure not consistent. 

The section shown in Fig. tla and b runs 
roughly in a north-east direction from the Con- 
nenta] Shelf up Spencer Gulf (C.S.1.R.0. 
Aust. L967a}. The gradient of the lines adja- 
cent to the rapidly deepening region outside 
Spencer Gulf under Station 64 again indicates 
a current inwards across the plane of the sec- 
tion, probably angled up slope so that the water 
inundates the more level sea bottom of Spencer 
Gulf, ‘This is again in agreement with the pre- 
viously deduced direction of flow into Spencer 
Gulf which is NNW and thus crosses the Jine 
of the section, 

Fig, 12 shows that the direction of the circu- 
lation in the Southern Ocean is consistent with 
the flow in the region outside the Gulf as de- 
duced from Figs 10 and 11, 


Conclusion 
The various data presented were mainly col- 
lected in the months of February to May, 
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; 11. Vertical sections along Spencer Gulf to 
the Continental Shelf. Marked —— on 
Figure 1, data taken from CSIRO Aust, 
1967a. Sampling depths indicated by 
dots. 

(a) Temperature (°C), (b) cr. 


although during different years, so that it would 
appear that the general circulation deduced 
occurs at least seasonally. However, to gain a 
complete and detailed piciure of the circulation 
in the Gulf, it is obviously necessary to con- 
duct a series of thorough cruises at various 
times throughout the year. In addition, direct 
current measurements need to be made ta 
check the results from the model. Since the dis- 


1-946 


4,923 


Fig. 12. Dynamic topography relative to 2000 db., 


south of Australia (taken from Bye 
1971}. 


tributions of water properties in St Vincent 
Gulf have been found (Bullock 1974) to be 
similar to those in Spencer Gulf, it may be 
expected that the general circulation in St 
Vincent Gulf would also display a pattern 
basically similar to that in Spencer Gulf. 
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